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The existence of phase transitions from liquid to gas phases in asymmetric nuclear matter (ANM) 
is related with the instability regions which are limited by the spinodals. In this work we investigate 
the instabilities in ANM described within relativistic mean field hadron models, both with constant 
and density dependent couplings at zero and finite temperatures. In calculating the proton and 
neutron chemical potentials we have used an expansion in terms of Bessel functions that is convenient 
at low densities. The role of the isovector scalar 5-meson is also investigated in the framework of 
relativistic mean field models and density dependent hadronic models. It is shown that the main 
differences occur at finite temperature and large isospin asymmetry close to the boundary of the 
instability regions. 
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I. INTRODUCTION 



Recently there has been a big development in the description of nuclei and nuclear matter in terms of relativistic 

many body theory. In particular, the phenomenological models developed using the relativistic mean field theory 

•^ . describe well the ground state of both stable and unstable nuclei [l|, |2|. These same models, with conveniently 

^»^ ' adjusted parameters, are used to describe the properties of neutron stars and supernovae |3, |j]. Therefore, it is 

f—^ , important to test these models at finite temperature and different densities. 

\^ The discussion of the properties of asymmetric nuclear matter (ANM) systems, namely their instabilities and phase 

C^ transitions, are presently a topic of great insterest [a, la B la JlOt Il3- The instabilities present in ANM may 

■^ manifest themselves as an isospin distillation or fractionation [l2|. It has recently been discussed [3, Ii3 that in 

^p \ an ANM system not only the mechanical but also the chemical instabilities appear as an instability of the system 

(^ t against isoscalar fluctuations. Hence the spinodal instability is dominated by density fluctuations which lead to a 

"'t' liquid gas separation with restoration of the isospin symmetry in the dense phase. This is known as the fractionation 

'~^ , effect. Multifragmentation also takes place when the system enters the spinodal region through nucleation or through 

r^ • spinodal decomposition. The correlation between spinodal decomposition and negative heat capacity evidences the 

r^ fact that the spinodal decomposition is the dynamics underlying the liquid gas phase transition [l3|. In the afore 

• • • mentioned work, the authors point out the fact that systematic measurements including correlations between these 

signals and also the inclusion of the isospin degree of freedom are important sources of future experimental work. 

On the other hand it is expected that neutron stars have a solid crust formed by non uniform neutron rich matter 

$^ in /^-equilibrium above a liquid mantle. In the inner crust, nuclei which form a lattice to reduce the Coulomb energy, 

coexist with a gas of neutrons which have dripped out. The solid crust has an important role in explaining the sudden 

spin jumps known as glitches in neutron stars |l4| . Recently it has been proposed that there is a relationship between 

the neutron skin of heavy nuclei and the properties of neutron star crusts |j|, namely the thicker the neutron skin 

of a heavy nucleus the thinner the solid crust of a neutron star. Properties of the crust, namely its thickness and 

pressure at the crust-core interface, depend largely on the density dependence of the equation of state (EOS). It is 

of particular importance the transition density below which uniform neutron rich matter becomes unstable. In 4] it 

was also shown that the thicker the neutron skin the lower the transition density form a uniform to a non uniform 

neutron rich matter. 

Within the framework of relativistic models, the liquid gas phase transition in nuclear matter has been investigated 
at zero and finite temperatures for symmetric and asymmetric infinite systems |a, llj, Ha . Il7| . With the help of 
the Thomas - Fermi approximation, we have investigated droplet formation in the liquid gas phase transition in 
cold ligl and hot |l3 asymmetric nuclear matters using relativistic mean field (RMF) models P, |23| • In |23| we have 
considered a droplet immersed in a gas of evaporated particles, in such a way that they mimic a source of changing 
mass. As temperature increases particles evaporate, mainly neutrons, and the fraction of protons in the droplet 
increases leading to isospin fractionation. In 22i| we have shown that at finite temperature droplet properties within 
different parametrizations of the RMF model have different behaviors with temperature. 

The EOS of neutron rich matter is in particular sensitive to the density dependence of the symmetry energy. In 
U the effect of changing this dependence was studied through the inclusion of non linear a ~ p and lo — p couplings. 



On the other hand the authors of [23 claim that the isovector scalar meson S is of vital importance in finding the 
stability conditions of drip-line exotic nuclei because the structure of relativistic interactions with a scalar (attractive) 
and a vector (repulsive) potential, which balance each other, is also present in the isovector channel. The J-field has 
also important effects on the symmetry properties of the nuclear system. 

Standard RMF interactions have their limitations even for describing nuclei close to the stablility line. This is 
due to the fact that the isovector channel is poorly constrained by experimental data. An example is the systematic 
overestimate of the neutron skins ^4] . Some of these limitations are overcome by quantum hadrodynamical models 
with density dependent meson-nuclcon couplings (which we refer to as TW) 25, 26, 27], which have been used 
with success to describe both nuclear matter and finite nuclei. In these models the couplings are either taken from 
Dirac-Brueckner-Hartree-Fock calculations or are fitted to data of nuclear matter or finite nuclei. 

We may ask whether the recent improvements of the RMF models, both through the inclusion of density dependent 
meson-nucleon couplings and/or the S scalar-isovector meson, present different features at subnuclear densities of 
nuclear asymmetric matter which could have consequences for the properties of the inner crust of neutron stars or 
in multifragmentation or isospin fractionation reactions. The parametrizations of these models take generally into 
account saturation properties of nuclear matter and properties of stable nuclei. Extension of the model for very 
asymmetric nuclear matter or to finite temperatures may show different behaviors. The simplest test between the 
models is a comparison of the regions of uniform unstable matter. 

In the present work we direct out investigation to some of the topics we have briefly mentioned before. We study the 
liquid gas phase transition and, in particular, determine the instability regions occuring within the density dependent 
hadronic models |2^ [23 and the non linear relativistic mean field models. We also study models including the S 
meson. 

In order to determine the instability regions, we found convenient to generalize the calculation of the chemical 
potentials as prescribed in (29j to relativisitc models with density dependent effective masses. The expansion we 
present in our paper is accurate in the range of temperatures and densities discussed and allows us to perform explicit 
derivatives on the chemical potentials. 

This paper is organized as follows: in section 2 we present the whole formalism in which our calculations are 
based, namely the relativistic models, the chemical potential expansions and the thcrmodynamical conditions for the 
description of the instability regions; in section 3 a brief explanation regarding the introduction of the isovector-scalar 
meson is reported; in section 4 the results are presented and the conclusions are drawn. Finally we include an appendix 
with some important formulae not given in the main text. 

II. THE FORMALISM 

We start from the lagrangian density of the relativistic TW model |25l | 

C = ij 



7^ ( za^ - T^V^ _ Dif . foP ) _ (M - Ts 



tP 



+ i(9^(/.9''0 - m2</)2) - In^M""" + \mlV^V^ - -^B^, ■ B^'' + ^m% ■ 6^ (1) 

where fl^^ — d^jVn — d^V^ , B^^^ — d^h^ — d^h^ — rp(fe^ x fe^). The parameters of the model are: the nucleon 
mass M — 939 MeV, the masses of the mesons nis, rriv, mp, and the density dependent coupling constants F^, r„ 
and Fp, which are adjusted in order to reproduce some of the nuclear matter bulk properties, using the following 
parametrization: 

^i{p) = 'ri{Psat)fi{x), i^s,v (2) 

with 

where x — p/psat and 

^p{p) = ^piPsat) exp[-ap(x - 1)], (4) 



with the values of the parameters rrii, Ti, Oubi, Ci and di, i = s,v,p given in |25|. Other possibihties for these 
parameters are also found in the literature J33- Notice that in this model the non-linear terms are not present, in 
contrast with the usual non- linear Walecka model (NLWM). 

Within the Thomas-Fermi approximation, the thermodynamic potential is obtained. After it is minimized with 
respect to the meson fields, the following equations are obtained 

00 = —^Ps, (5) 



Vo - ^P, (6) 

mf, 



bo = ^P3, (7) 



with 



d^p 
P^Pp + Pn, P3^Pp-p„, Pi == 2 / —--^(/i+ -/.,_), i=p,n (8) 

Ps — Psp + Psn, Psi = 2 / \3,~j^^fi+ + •^^-)' 



M* =M -Ts 00, E* = Vp2+M*2 and /,± = 1/{1 + cxp[(^* =F i'i)/T]} , where the effective chemical potential is 

p 

u^ = fi^~ TyVo - Ta-^bfi - S^, Tp3 = 1, r„3 = -1 (9) 



with the rearrangement term given by 



j:^=—pVo + — p,-- — p.^o. (10) 



The energy density in the mean field approximation reads: 



-'^0 + ^r^o +^^%. (11) 



2 "" 2 " 2 
The pressure becomes 

p dT, 



p(r„r„,rp) = -L ^ y dp|^ (/.+ + /._)- ^^^ |^i + 2 



Ts dp 



For two of the usual NLWM parametrizations, namely NL3 ^1] and TMl '3^, the above equations read: 

£iVL3 = ^(5s,5.,5p) + ^ + ^, PNL3^P{gs.g..gp)-'^-^ (13) 



and 



•tTMi -t(gs,gv,gp) + — + -^^ H g — , -^ta/i - P[gs,gv,gp) g 24" "' 24~' ^ '^ 



where the meson- nucleon couphng constants, Qs, Qv, and gp substitute rsjr^,, and Tp. They are not density dependent 
and consequently all derivative terms in the pressure cancel out and k, A and i^ are the self-coupling constants 
multiplying the non-linear terms. 

In order to study the instability region at low densities it is convenient to invert eq.® and obtain the effective 
chemical potential that appear inside the distribution functions. We have followed the prescription given in 6. 29|. 
where just the particle distribution function is considered. The effective chemical potential reads: 




where 



7 is the spin multiplicity, /3 = 1/T and 



E^d 



jQiPM* 



d^'p ^_pE' _ M 



■ p,n (15) 
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with Kn the modified Bessel function. The coefficients bi are defined in terms of the ratios 

_ Q{n(3M*) _ K2{nf3M*) 
"~ QiPM*) ~ nK2{l3M*y 

and have been calculated in [23. We list the first three 

5i = 02, 62 = 03 - a2/2, 63 = 04-0203 + 02/3, 

with 02 = 5'2(/3M*), 03 = 2S'|(/3M*)-53(/3M*), a^ ^ bSl{l3 M*) ~ b S^ifi M*)S2{P M*) + S^ifi M*). The chemical 
potentials /Zj , i = p,n are obtained from equation Q and they read 



= r„ K, + ns^bo + S^^ + - ln(770 + > ^ birji (17) 




where the rearrangement term is given in equation (|10|l . 

We have checked the range of applicability of the expansion given in eg. ljlSII by comparing it with the exact results. 
We have concluded that, for symmetric matter, the expansion works very well for temperatures T > 7 MeV and 
subsaturation densities, p < po. This agreement improves for higher temperatures within the range of temperatures 
involved in the present work, where the antiparticles do not play a crucial role. For ANM, T =7 MeV, the agreement 
is still good for the proton chemical potential. The neutron chemical potential is well reproduced only for p < O.Spo as 
can be seen in fig. ^ where we have plotted the effective chemical potentials for protons and neutrons for asymmetric 
matter with the proton fraction j/p = 0.1, where yp = —. This result is still adequate for the range of densities we 
have considered in the present work 

The stability conditions for asymmetric nuclear matter, keeping volume and temperature constant, are obtained 
from the free energy density T , imposing that this function is a convex function of the densities pp and p„, i.e. the 
symmetric matrix with the elements [5|, [3, llOJ 






dpidpj 



is positive. This is equivalent to imposing |3S 



P^ > 0, (18) 

dpp 

P^ > 0, (19) 

d{Pp, Pn) 



where we have used ^i = 4^ 
respectively, in the form 



and 



In terms of the proton fraction the conditions H18|) and (|19|l can be rewritten, 



'Pj^ti 



dP 

dp 



>0 



,yp 



dP 



, (^] > 0. 



(20) 



It has recently been argued *3, 'TOj that in ANM the spinodal instabilities cannot be separately classified as mechanical 
or chemical instabilities. In fact, the two conditions that give rise to the instability of the system are coupled so that 
it appears as a mixture of baryon density and concentration fluctuations. Therefore, we define the stability region as 
determined by iPHjI . 

In order to calculate the boundaries of the spinodal instability regions we use the Gibbs-Duhem relation, at a fixed 
temperature and isospin asymmetry. 
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dp 



with 5 = —pzjp = 1 — 2yp, and 



dp, 
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dp 
'dS 
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dP 
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T,p 



where 



dP 

'dS 



T,p 



n ^x\^^'^ ^f^ x\^^p 



El dpi 



t^p,n 



(21) 



(22) 



(23) 



With the expressions given in this section the spinodal regions can be obtained for both different temperatures and 
different parametrizations. 



III. INCLUDING ISOVECTOR-SCALAR MESONS 



To investigate the influence of the ^-meson in the stability conditions we have included in the NLWM the isovector- 
scalar meson terms l23l: 



£ = '4^ 



[u 



9p, 



^^{td''-g.V^~^f-b^)-{M-gs<j>-gsT-S) 



v- 



+ 1 (5^09'^0 - m202 _ 1^^3 _ _L^^4\ _ l^^^f^M- + ImlV^V^ 



lB^,.B'^- + lm%.b^ 



-{dJd^5~mlP), 



(24) 



where gg and m^ are respectively the coupling constant of the 5 meson with the nucleons and its mass. Self-interacting 
terms for the u-meson are also included, k, and A denoting the corresponding coupling constants. The set of constants 

is defined by g^ = ^/.mf, i = s,v,5,^^ J Ipfn% rris = 550 MeV, m^, = 783 MeV, rup = 763 MeV, ms = 980 MeV, 

f, = 10.33 fm2, /„ = 5.42 fm^, fp = 3.15 fm^, fs = 2.5 fm^, k = 0.0665^ and A = -6 x 0.0048^4 ^ and we call 
it NL(5. For reference, in Tabled we show the properties of nuclear matter reproduced by the models we discuss in 
the present work. From the minimization of the thermodynamic potential obtained in a Thomas- Fermi approach, the 
equation of motion for this field becomes 



X 9s 

03 = -^Ps3, 



(25) 



with ps3 = Psp — psn ■ The energy density and the pressure are also affected by the presence of the new meson. The 
term +1/2 to^ (5| should be added to the energy density and the — l/2m^(5| should be added to the expression of the 
pressure, both given in eq. ^\?>\ . The effective masses for protons and neutrons acquire different values, namely. 



M* = M - gs (l)a - Tagsh 



p,n. 



For completeness, we have also included the (5-meson in a model where the p and 6 couplings are density dependent, 
as done in 2^1, where it is called density dependent hadronic model (DDHpd) . For this purpose, we have considered 
the density dependence of the p and (5-nucleon vertices given in fig. 1 of ref. zS\ which have been extracted from DBHF 
calculations of ref. |34|. In this case the coupling constants gs, gv, 5p and gs used in equation H24|l should be replaced 
by Fg, F^,, Tp and F^ and the non- linear scalar terms do not appear. For F^ and F,y we take the parametrizations 
given in Eqs.Q and (@J. For Tp and F^, we propose the following paramctrization. 



Ti = T{psat) Mx), X = p/ps 



with 



f{x) = aiex:i>[-bi{x - 1)] - Ci{x - di), 



P, 



and the parameters ai,bi,Ci and di are defined in Table ITU This paramctrization reproduces the curves given in fig. 
1 of [23 and is also displayed in fig. |21of the present work. We point out that in the present work there is a factor 2 
difference in the definition of p- meson coupling constant. We take for the effective chemical potential Eq. 0, with 
the rearrangement term given by 



and for the pressure 



dp 



pVo 



dTp 60 dTs 



dp 



P3y 



dp 



Ps4>Q 



dTs 
dp 



PS3S3 



P(F„F„Fp,F5)-P(F„F„Fp)--^53^ 1 



, P dTs 
'Ts dp 



IV. RESULTS AND CONCLUSIONS 

A quantity of interest in ANM is the nuclear bulk symmetry energy discussed in |23j . This quantity is important 
in studies involving neutron stars and radioactive nuclei. The behavior of the symmetry energy at densities larger 
than nuclear saturation density is still not well established. In general, relativistic and non-relativistic models give 
different predictions for the symmetry energy. 

It is usually defined as 



_ 1 d^£/p 



(26) 



s=o 



which, for the models not including the (5-meson, can be analytically rewritten as 



^svni — 






Tl 
8^' 



(27) 



where 

with kp = (I.Stt^p)^/'^. If the S meson is included the symmetry energy reads [2 



1/3 



k^ 



F2 

P 



Tl /M* 



\E*f 



(28) 



In fig|3|we show the symmetry energy for the different models used in this work, calculated at T = and j/p = 0.5. 
We can see that at subsaturation densities where the instability regions occur (0 up to 0.1 fm"'^), NL3 and TMl 



parametrizations give very similar behaviours. The NLS parametrization describes saturation at 0.16 fm"'' and this 
fact affects the symnietry energy at low densities, in particular it is the model with the lowest symmetry energy in this 
range of densities (see Table QJ. The TW model presents a different behaviour, after a faster increase at low densities, 
the symmetry energy increases much slower than most of the other models at larger densities. A very low value of 
the symmetry energy was obtained with the DDHpd parametrization. It amounts to 25 MeV, in contrast with the 
value mentioned in 28] (33.4 MeV). Nevertheless, in the same reference |23, the authors show a curve in fig. 2, which 
confirms our number. The density dependent hadronic models show a softer symmetry energy at the densities shown. 
At higher densities the models containing the (5-meson are expected to get a harder behaviour due to relativistic effects, 
namely the S contribution goes to zero and only the repulsive contribution from the p-meson remains |2^. We point 
out that in present work we are only testing the low density region of the symmetry energy, namely our discussion is 
concerned with the behaviour of nuclear matter at densities below O.lpQ, where po is the saturation density. In this 



9£a 



and curvature 

P=PO 



region the quantities which better distinguish the different models are the slope L = Spo I — g- 

Ksyrn = 9po I Q^T^ ) of the symmetry energy 35] given in Table ^ The TW model is the only one with a 

P=PO 

negative curvature at saturation density. In fact, the DDHp^ also has a negative curvature at densites below 0.13 
fm~'^ but at saturation the curvature has already a positive value. 

In fig0]we plot the spinodals for the NL3 parameter set and different temperatures, varying from T=0 up to T=14 
MeV, already close to the critical temperature, Tnls.c ^ 15 MeV. All the other models present a similar behaviour, 
namely the larger the temperature the smaller the instability region. These same conclusions have already been 
discussed in [g for Skyrme type phenomenologial EOS for ANM. After a critical temperature the liquid-gas phase 
transition is a smooth transition. A consequence of this behavior is the greater stability against density flucutations 
of neutron-rich systems. Also, in intermediate energy collisions fragmentation or fractionation occurs later after the 
system has cooled down. 

In fig. |31we display the spinodals for different parameter sets obtained with T —0, 10, 14 MeV. For T = MeV 
the differences are not significant, occuring at the higher density branch, ^ 0.1 fm~'^. In particular, in the NL(5 
parametrization the boundary lies at a larger density than in the other models for S < 0.5. The spinodals for the 
density dependent hadronic models are decreasing less with the asymmetry parameter S than the others. At finite 
temperature the differences are larger. These differences occur again in the larger density branch, and it is again 
the NL(5 which presents a boundary at larger densities in the small asymmetry region followed by the TW and the 
DDUpS models. In the large asymmetry region we are already testing both the critical temperatures and asymmetries 
of the models and the differences between the models are larger. For the DDHpJ model, as temperature increases, the 
instability boundaries extend to higher asymmetries as compared with all the other models. This is possibly due to 
a lower symmetry energy. The opposite occurs with the NL3 parametrization which has the larger symmetry energy. 
We conclude that the information obtained from the spinodal decomposition sensitive to the underlying model comes 
from the phase space close to the spinodal boundary at large isospin asymmetry. 

In conclusion, we have compared the spinodal boundary as a function of density and isospin asymmetry, at sev- 
eral temperatures, for different relativistic models. We have considered both quantum hadrodynamical approaches 
with non-linear terms (NL3, TMl, NLJ) and density dependent hadronic models with density dependent coupling 
parameters (TW, T)DIlpS). The largest differences between the models occur at finite temperature and are more 
clearly shown in the high isospin asymmetry region. Only physic al q uantities that explore this region in phase space, 
as for example the neutron-proton differential fiow suggested in |36|. will bring possible information about the EOS 
of ANM. Another physical system that could be sensitive to present results is a protoneutron star in the process of 
stellar colapse. The physical structure of the matter, namely a smaller or larger non-uniform nuclear matter region, 
may affect properties such as the neutrino scattering processes |37|. 

Another point to be investigated is the role of the Coulomb interaction and finite size effects in the instabilities 
of ANM. We have already shown that the electromagnetic force cannot be neglected in nucleation processes, where 
droplets of lower asymmetry are formed in a very asymmetric gas background either at T = or finite temperatures 
[isL [13 • The authors of ref. Q have included the Coulomb interaction and surface tension in their calculation of 
the binodals in a simplified way and they have shown that the minimum pressure for a given temperature do not 
occur for symmetric nuclear matter and that surface effects lower the binodal pressure. It has also been shown that 
the Coulomb interaction affects the growth of instabilities |33] . A natural extension of the present work would be 
the inclusion of the Coulomb interaction in the same spirit as done in 9] and consequent construction of the binodal 
sections and spinodal regions for different temperatures. We expect, however, that the overall effect of the inclusion 
of the Coulomb interaction will be the same for all the models and the main conclusions of the present work will not 
change. 



V. APPENDIX: SOME FORMULAE 

In what follows T is considered fixed. We give the main expressions needed to obtain the conditions which define 
the boundaries of the instability regions. We first consider the derivatives of the chemical potential with respect to p, 
for a given S. For this purpose we need the derivative of the effective chemical potential that reads: 



di^i \ dvi drji 1 Y^ ^ dhk dM* 



where 



di^t _ 1 

% ~ 13 



^+Efc^^^'" 



and 



P 2 

So, in the NLWM model, one has 



dru^ru_ f3^ ^^^p^ [4X2(x) - x {Ki{x) + K^{x)) 
dp 



x^Kl{x) 



dM* 



x^PM* 



djM_ 
9p / s 



dui 



dp / s 
Instead, for the TW model, the previous equation reads 



Tj3 






dp J s \dp ) s dp\ml J '^9pi4m2 . 



dp ' 



with 



S^^ = 



r„ dV, 2 , Tp dTp 2 .2 rnl dT, 



ml dp Am? dp 



^'^ Tf^(^-^'*) ' 



M - M* = — p, 
m. 



To evaluate the expressions involved in H22|l we obtain the derivative of the chemical potential with respect to 5 for 
a fixed p 

dvi \ _ dvi d^^ ^ 1 -^ ^ dbk dM* 



dS J ^ dri, d5 (3^ " dM* dS ' 

^ k 



where 



T3iP P"^ TT^Pi 



AK2{x)-x{Ki{x)+K^{x)) 



drg 

dS ^ 2-fQ{x) 2 [ x^K^ix) 

Both the derivatives ^|^ and ^^|- are calculated numerically. 



dM* 
"df 



The expressions for ( -^ j become 



dpi 
dS 



dvi 
'dS 



^'' (2I; ' ' 



in the NLWM, and 



'^^ / p \ '^^ J p V '^''^p J 2 ml dp 86 



2 



d,.\ _ fd^A rJllL.r,+ Il^9JlL,H,9^^ 



in the TW model. 

Some derivatives of the coefficients bi are also listed: 



dbi dS2 

P- 



dM* dx 



db2 pdh,^/^^^dS, dS, 



dM* dx \ dx dx 

with 

dSn{x) 1 



dx 2niir|(x) 



[K2{nx) {Ki{x) + Kaix)) -nK^ix) {Ki{nx) + Ksinx))] 
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TABLE I: Nuclear matter properties. 





NL3 f31] 


TMl [32] 


TW [25] 


NL5 [23] 


DDHp5 [22] 


B/A (MeV) 


16.3 


16.3 


16.3 


16.0 


16.3 


Po (fm-^) 


0.148 


0.145 


0.153 


0.160 


0.153 


K (MeV) 


272 


281 


240 


240 


240 


Ssyr.. (MeV) 


37.4 


36.9 


32.0 


30.5 


25.1 


W/M 


0.60 


0.63 


0.56 


0.60 


0.56 


L (MeV) 


123. 


117. 


55. 


101. 


44. 


Ksy,^ (MeV) 


108. 


36 


-124. 


112. 


45. 



TABLE IL Parameters of the DDHpS model 



i Fi 


fli 


bi 


Ci 


d. 


p 5.8635 


0.095268 


2.171 


0.05336 


17.8431 


5 7.58963 


0.01984 


3.4732 


-0.0908 


-9.811 
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FIG. 1: Chemical potential curves for T — 7 MeV, NL3 parameter set and j/p = 0.1. The curves for the exact proton chemical 
potential and its expansion are coincident (short dashed line). The straight long dashed line represents the exact neutron 
chemical potential and the bending one (dotted line) its expansion. The expansions are taken up to S*** order. 
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FIG. 2: Density de pen dent couplings as obtained from the proposed parametrization given in the text and which reproduce 
the curves given in J2^. 
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FIG. 3: Symmetry energy results for the NL3, TMI, TW, NLJ and DDHp5 parameter sets. 
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FIG. 4: Spinodal regions obtained with the NL3 parametrization and different temperatures. The instability regions lie inside 
the curves. The curves are drawn for T=0 (solid line), 8 (dotted line), 10 (short-dashed line), 12 (dash-dotted line) and 14 
(long dashed line) MeV. 
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FIG. 5: Spinodal regions for different parameter sets and T = (upper panel), T = 10 MeV (middle panel) and T = 14 MeV 
(lower panel). 



